extern doublereal enorm_(integer *n, doublereal *x);
